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Abstract: Extragalactic surveys in the emission line of [O III] A5007 have provided us with the 
absolute line strengths of large, homogeneous sets of planetary nebulae. These data have been used 
to address a host of problems, from the measurement of the extragalactic distance scale, to the study 
of stellar populations. I review our current understanding of the [O III] planetary nebula luminosity 
function (PNLF), and discuss some of the physical processes that effect its structure. 1 also describe 
the features of the Ha PNLF, a function that, upon first glance, looks similar to the [O III] PNLF, but 
which includes a very different set of objects. Finally, 1 discuss recent measurements of a, the number 
of PNe found in a stellar population, normalized to that population's bolometric luminosity. I show 
that, contrary to expectations, the values of a found in actively star-forming spirals is essentially the 
same as those measured in late- type elliptical and lenticular systems. I discuss how this result sheds 
light on the physics of the planetary nebula phenomenon. 
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1 Introduction 

Photometric surveys for extragalactic planetary nebu- 
lae (PNe) are useful for a host of reasons: they produce 
test particles for studies of galactic kinematics and 
chemical evolution, they probe facets of stellar evo- 
lution that are unobservable with conventional tech- 
niques, and they yield unique insights into the proper- 
ties of Galactic PNe. PN photometry even allows us to 
measure the extragalactic distance scale to a precision 
that few other methods can reach. These uses come 
about because, unlike Galactic samples of PNe, the 
planetary nebulae of other galaxies are all at the same 
distance and all have the same amount of foreground 
reddening. Thus, we have access to the objects' ab- 
solute line-strengths, and, in particular, the planetary 
nebula luminosity function (PNLF). 

The hist ory of the PNLF date s back to the early 
1960's, whenlHe nize fc Westerlundl (|1963D modeled the 
PNe of the Magellanic Clouds as a set of spherical gas 
clouds surrounding non-evolving central stars. Under 
this assumption, the expected Balmer-line flux from a 
PN should go as 



F oc NhN^V oc oc t~ 



(1) 



which, when written in terms of magnitudes, yields 

t oc IQ"^'-^ cx 6°-^"^'^^ (2) 

Since the number of objects observed between mag- 
nitudes M and M -\- dM is proportional to the time 
spent between those magnitudes, the implied luminos- 
ity function of Magellanic Cloud PNe is then 



N(M) oc — oc e"-^"™ 
^ ' dM 



(3) 



Obviously, the above argument grossly oversimplifies 
the physics of the planetary nebula phenomenon. Most 



PNe are not spherically symmetric, and their early evo- 
lution is governed more by the hydrodynamics of a fast 
wind/sup erwind inter action than by simple photoion- 
ization and expansion j|Kwokll982l : lBalick fc Frankl2002l : 
ISchonberner et"al ] |2007D . As a result, the ionized mass 
of a young PN should increase with time and produce 
a distribution of PN luminosities that is flatter than 
that implied by equation ((3)1 . Furthermore, as a PN 
evolves, its opacity to ionizing photons and its ioniza- 
tion structure may change as well, adding dips and 
wiggles to the expecte d form of the PNLF. Yet for all 
this, the Jacoby (1980) observations of faint PNe in the 
Large M agellanic Clou d seem ed to confirm the simple 
iHenize fc Westerlundl (jl963l l law. Although the LMC 
observations were conducted at 5007 A, rather than 
Ha (since PNe are easier to identify in the [O III] line), 
the agreement between the observations and the sim- 
ple theory suggested that much of the complex physics 
associated with planetary nebulae was not crucial to 
understanding their faint, late stages of evolution. 

The late 1980's saw the introduction of the [O III] 
A5007 PNLF as an extragalactic standard candle. Sur- 
veys of more than two dozen spiral and elliptical galax- 
ies demonstrated that, while the faint-end of the PNLF 
may be fit by a Henizc fc^estcrlund (1963) exponen- 
tial, the bright-end of th e func tion truncates abruptly. 
This led ICiarduUo et~al] (jl989i l to propose fitting the 
PNLF using a cutoff exponential. 



N{M) 



1 - e 



(4) 



where M* = —4.48 is the absolute luminosity of the 
brightest possible planetary in a magnitude system de- 
fined by 



M5007 = -2.51ogF5oo7 - 13.74 



(5) 



Remarkably, Af* appeared to be virtually independent 
of stellar population: internal tests within galaxies 
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Figure 1: Values of Af* derived using the distances of 
13 galaxies with Cepheid measurements. No correction 
for metallicity has been applied. The abscissa shows 
the galaxies' oxyge n abundances, as est imated from 
their H II regions l|Ferrarese et al.l 120001 ). The curve 
show s the theoretical res ponse of the nebula to metal- 
licity [DoEitZ^eF^l] (jl992l l. modified to reflect the pres- 
ence of metal-poor stars in metal-rich galaxies. The 
consistency of the measurements is excellent, and the 
scatter agrees with the internal errors of the methods. 




(e.g., ICiarduUd l2003l : iHui et all 1 19931 ) internal tests 
withi n galaxy gr oups and clusters (e.g..'Ciard ullo et al.l 
Il989l : [jacobv etall 1 19901 : [Ciarduho et al.. ,19"9ll ). and 
external tests in galaxies with distances from Ce pheid 



Ionization Magnitude 

Figure 2: The planetary nebula luminosity func- 
tions expected from nebulae which reprocess their 
central star's ionizing radiation into [O III] A5007 
flux at 10% efficiency. The histograms are based on 
the post-AGB hydrogen-bu rning evolutionary tracks of 
IVassiUadis fc WoodI (fl994l ). Note the deficit of objects 
at intermediate magnitudes: this is due to the rapid 
cooling that occurs immediately after shell-burning 
ceases. The figure also suggests that for the highest- 
mass cores, the conversion efficiency of ionizing flux to 
[O III] A5007 radiation in no more than a few percent. 



and s urface brightness fluctuation measurements (jCiardullo 
|2002| ) yielded values of M* that were consistent to bet- 
ter than ~ 10%. As Figure [1] illustrates, the only sys- 
tematic trend was a slight fading of M* at low metal- 
licity - a fading that was entirely consistent with the 
expected response of a set of nebulae to lower oxygen 
abundance (C iarduU o fc Jacobvll9"9^ : iDopita et al.ll993 : 
ICiardullo et al.ii200a ). 



2 The Faint-End Shape 

The constancy of M* means that the bright-end cutoff 
of the [O III] A5007 PNLF contains little or no infor- 
mation about the underlying stellar population. How- 
ever, the same is not necessarily true for the function's 
faint-end shape. As described above, the physics of FN 
evolution is extremely complicated, and a proper treat- 
ment of the pro blem requires careful hy drodynamical 
modeling (e.g.. ISchonberner et al.ll2007l ). Yet at its 
most fundamental level, FN evolution is controlled by 
two timescales: one tied to central star evolution, and 
the other to nebular expansion. If most FNe are pow- 
ered by low-mass, slowly-evolving central stars, then 
nebular expansion will dominate and the result will be 
a luminosity function that monotonically increases to- 
wards fainter magnitudes. Conversely, if the majority 
of central stars are high-mass objects, then their rapid 
stellar evolution will render nebular expansion unim- 
portant. In this case, the observed FNLF will primar- 
ily reflect the evolution of the central star's ionizing 
flux, and there will be a deflcit of FNe with luminosi- 



et al.l 



ties near the top of the white-dwarf cooling sequence. 

Figure [2] illustrates this effec t by trans lating the 
post- AGB evolutionary tracks of IVassiliadi s & Woq^ 
(jl994l ) into a set of mock luminosity functions, un- 
der the simplistic assumption that planetary nebu- 
lae reprocess ionizing flux into [O III] A5007 radiation 
with 10% efficiency. Obviously, many processes can re- 
duce this factor, and real FNe may evolve quite differ- 
ently. Nevertheless, more sophisticated analyses bare 
out the result of this straightforward analysis: popula- 
tions which are dominated by high core mass FNe must 
have a dip in their luminosity function (jMendez et al.l 
120081 ). 

The stars of M31's bulge have ages of ^ 6 Gyr 
(jTrager et al.l [2000) . C onsequently, from the initial 
mass- final mass relation (|Kalirai et al ]|2008l ). we would 
expect the region to have FNe with low-mass, slowly 
evolving central stars. This is i ndeed the case: the 
deep [O III] luminosity function of lCiardullo et al.l (|2002l ) 
shows that the population' s FNLF is well-fit by the 
iHenize fc Westerlundl (( l963h law for non-evolving cen- 
tral stars. Conversely, the star- forming populations 
of th e Small Magellanic Cloud llJacobv fc De Marcol 
|2002| ) and M33 (|Ciardullo et aLll2004l ) should have manv 
high-mass cores that evolve rapidly across the HR di- 
agram. Again, this is what is seen, as the luminosity 
functions of both systems have a clear deficit of inter- 
mediate luminosity objects. 

How fast is this transition? In theory, we can 
constrain this quantity by examining the FNLF in 
a galaxy where star-formation has ceased in the re- 



www. publish . csiro. au /journals /pasa 



3 



C! 3 



E 10 



M31 Bulge 
(U-V)„ = 1.54 



NGC 5128 
(U-V)„= 1.29 



M33 

(U-V)„ = 0.31^ (k. 



V SMC ^ 



(U-V)„ = 0.18 

_l I I 



-4 -2 0-4-2 

Absolute A,5007 Magnitude 

Figure 3: Deep [O III] PNLFs for the bulge of M31, 
the envelope of NGC 5128, the disk of M33, and the 
Small Magellanic Cloud. Open circles denote bins ef- 
fected by incompleteness; t he dotted line shows the 
iHenize &: Westerlundl (119631 ) exponential law. Only 
the strongly star-forming systems have a deficit of PNe 
at intermediate magnitudes. 



cent past. Unfortunately, such objects are difficult 
to identify. One possible location is the envelope of 
NGC 5128: this interacting elliptical has a blue color 
that is reminiscent of the inter-arm region of a spi- 
ral galaxy, and still contains some e vidence of recent 
star- formation (|Reikuba et al.|[200l '). As Figure [3] il- 
lustrates, a deep planetary nebula survey of the re- 
gion reveals a PNLF that is identical to that found 
for M31's bulge; there is no hint of the mid-magnitude 
deficit present in the SMC and M33. This suggests 
that the transition from a central-star driven lumi- 
nosity function to a nebula-driven luminosity function 
must occur rather quickly: unless the PN population 
is d ominated by objects f r om ve ry young progenitors, 
the lHenize fc Westerlundl l| 19631 ) law is the better pre- 
dictor of the PNLF's shape. 



3 The PNLF in Ha 

Since the brightest PNe emit predominantly in [O III] 
A5007, it is natural that PN surveys in distant galax- 
ies focus on this wavelength. But other emission lines 
contain information as well, and for many Galactic sur- 
veys, it is the Ha emission that is most-easily studied. 
Thus, we should consider the probative value of this 
brightest Balmer line as well. 

Figure |4] displays the Ha luminosity function for 
three Local Group galaxies, using an Ha magnitude 
defined via 



Mhc = -2.51og_FHc - 13.67 



(6) 



Like the [O III] PNLF, the Ha luminosity function dis- 
plays a cutofi' that is insensitive to stellar population. 
The planetaries produced by the old stellar popula- 
tion of M31's bulge have roughly the same maximum 
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Figure 4: The Ha PNLFs for M31's bulge, M33's disk 
and the LMC. Though the shape of the cutoff differs 
from galaxy to galaxy, the brightest PN in each system 
has about the same absolute brightness. The turnover 
at faint magnitudes is partially due to incompleteness, 
as we are using [O III] selected samples. 



Ha fiux as the PNe created in the star-forming en- 
vironments of M33 and the Large Magellanic Cloud. 
Although the cutoff in the Ha luminosity function is a 
factor of ~ 2.5 fainter than that for [O III] A5007, the 
data do suggest that Ha measurements of PN ensem- 
bles can be used to estimate distance. 

This is not to say that all the PNe that are bright 
in Ha are also bright in [O III]. In fact, when com- 
pared against the most luminous objects in the [O III] 
A5007 PNLF, Ha-bright objects are surprisingly inho- 
mogeneous. This can be seen in Figure O which gives 
the [O III] to Ha-|-[N II] line ratios of PNe in M31's 
bulge, M33's disk, and throughout the Large Magel- 
lanic Cloud. When plotted against [O III] absolute 
magnitude, the data show excellent consistency at the 
bright end, as the line-ratios, R, populate a wedge in 
emission-line space defined by 



4 > i? > 3.14 



(7) 



where L* = 2 .4 x 10 ergs s in t he monochromatic 
5007 A line l|Herrmann et al.ll2008l ). All PNe in the 
top magnitude of the luminosity function have [O III] 
brighter than Ha by a considerable amount. Con- 
versely, PNe in the top magnitude of the Ha PNLF 
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Figure 5: The ratio of [O III] A5007 to Ha+[N II] 
observed for PNe in M31's bulge, M33's disk, and the 
Large Magellanic Cloud. When plotted against [O III] 
brightness, the samples appear rather homogeneous, 
especially at the bright end. When plotted against Ha 
brightness, this consistency disappears. 



have line ratios that span the entire range observed, 
from 0.1 < 7i < 4. This scatter, which ha s been repro- 
duced in the hydrodynamic simulations of lMendez et al.l 
(|2008l ). clearly demonstrates that more than one type 
of object is contributing to the bright-end of this func- 
tion. 

Figure [5] also illustrates a danger about drawing 
conclusions from samples of objects selected via a sin- 
gle emission-line. PNe in distant galaxies are typi- 
cally identified from their flux at [O III] A5007; conse- 
quently, objects that are faint in [O III] but bright in 
Ha may not be detected. The result can be a censored 
dataset. So, while the top panel of Figure [S] is an ac- 
curate representation of the distribution of objects in 
emission-line space, the bottom panel of the figure is 
seriously incomplete at the faint end. The same is true 
for the luminosity functions of Figure |4l although the 
deficits of PNe at faint magnitudes may be real, sam- 
ple incompleteness is also playing an important role. 
In the Milky Way, where large sa mples of PNe are se- 
lected via their brightness in Ha (|Parker et al.ll2006l : 
iGonzalez-SoIares et all I2008D . this bias will be seen 
when objects are plotted against their [O III] A5007 
magnitude. 
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Figure 6: The observed ratio of [O HI] A5007 to 
Ha-f [N II] plotted against [O HI] A5007 absolute mag- 
nitude for those bright PNe in the LMC and M31 with 
published spectrophotometry. The arrows display the 
effects of internal extinction, as estimated from the 
objects' Balmer decrement. The location of the PNLF 
cutoff, M* , is identified via the dotted line. Note that, 
although no PN is ever observed to have an [O III] 
A5007 magnitude brighter than M* , many have in- 
trinsic luminosities much brighter than this value. 

4 The PNLF and Distances in 
the Milky Way 

Figures [31 |H and [S] suggest a method of constrain- 
ing the distances to PNe within the Milky Way. As 
the figures illustrate, PNe populate a specific region of 
emission-line space: there is a hard upper limit to a 
PN's luminosity, both in [O III] A5007 and Ha. More- 
over, objects near the bright-end cutoff of the [O III] lu- 
minosity function are all high-excitation objects. Thus 
the line strengths of a Galactic planetary nebula im- 
mediately provide an upper limit to distance. For ex- 
ample, a PN which is twice as bright in Ha than in 
[O III] cannot be any brighter than M5007 ~ — 2. 

Unfortunately, such relations are of limited use. 
Like all objects, PNe are affected by extinction along 
the line-of-sight. This extinction has two components: 
one from the foreground dust of the Milky Way Galaxy, 
and the other from circumstellar material. The effects 
of the former have to be removed, since the dust is to- 
tally unrelated to the planetary. Conversely, the latter 
extinction must remain untouched, as it is an integral 
part of the planetary nebula phenomenon. Indeed, as 
Figure |6] illustrates, circumstellar extinction helps de- 
fine the PNLF cutoff: without this dust, many of the 
PNe observed in the Large Magellanic Cloud and M31 
would be substantially brighter than AI* . Thus, if one 
wants to use PN line strengths as a distance indicator, 
one has to somehow separate the Galactic and circum- 
stellar components of extinction. 

5 The PNLF's NormaUzation 

An under-appreciated component of the PNLF is its 
normalization. According to the theory of stellar en- 
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ergy generation, all non-star forming stellar popula- 
tions should create B ~ 2 x 10~^^ PNe yr~^ Lq^, re- 
gardless of the s ystem's age, metall icity, or initial-mass 
function (Renzi ni &: Buzzonilll98(^ . Therefore, if one 
has an estimate for the mean PN lifetime, that age can 
immediately be translated into an expected number of 
observable planetaries. Any departure from this num- 
ber provides a constraint on the timescale of the PN 
phase and/or the importance of alternative channels 
of stellar evolution. 

To perform this experiment, one needs to estimate 
a, the number of PNe divided by the total bolometric 
luminosity of the underlying stellar population. This 
presents two problems. First, one has to translate 
broadband measurements of a galaxy's flux into to- 
tal bolometric luminosity. This is not as difficult as 
it sounds: the PNLF provides the galaxy's distance, 
and the population's bolometric correction can be as- 
sumed to be B.C. ~ —0.85. inde pendent of its age 
or metallicity (|Buzzoni et al.l [20061 ). Second, one can 
never observe the entire PN population of a distant 
galaxy: typically, one only counts objects in the top 
~ 1 mag of the luminosity function. Thus, extragalac- 
tic PN surveys usually quote ao.s or 02.5, i.e., the nor- 
malized number of PNe in the top 0.5 or 2.5 mag of 
the luminosity function. To obtain the total PN pop- 
ulation, this value has to be extrapolated over a range 
of at least 5 or 6 magnitudes. 

For examp le, in their 19 89 survey of the Leo I 
Cloud, Ciardul loeraLl ||1989D found values of 02.5 ~ 
40 X 10~^ PNe Lg^ for three different gala xies. They 
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Figure 7: The left hand panel displays values of ao.5 
for late-type stellar populations; the right hand panel 
plots a measurements performed in the inter-arm re- 
gions of spiral disks. Although there are trends in the 
data — large, red ellipticals are less efficient at pro- 
ducing [O IIl]-bright PNe than smaller, blue galaxies 
— the most striking feature is the similarity in the 
numbers. All stellar populations have about 1 bright 
PN for every 5 x 10* solar (bolometric) luminosities. 



On the other hand, estimates of the total PN pop- 
ulation require an extrapolation. For simplicity, this 
is usually performed by integrating a a fit to equa- 
tion 0. However, as Figure [3] shows, star-forming 
systems have a deficit of objects at intermediate mag- 
nitudes; if one integrates the SMC's PNLF found by 
[jacobv fc De Marcol l|2002l '). then one obtains a total 



then e xtrapolated this value to os.o using the lHenize fc WestArMuR^ Pulation that is ~ 1.6 times less than that de- 
(Il963l ) law, and assumed a timescale for the PN phe- 
nomenon of ~ 25,000 vr iPottaschI (119841 ). The result 
was a stellar death rate that was surprisingly close to 
the theoretical value for single stars. 

Unfo rtunately, the true situa tion is a bit more com- 
plicated. ICiarduUo et al.l (|2005t ) have presented values 
of ao.5 for ~ 20 elliptical and lenticular galaxies. As 
the left-hand panel of Figure [7] demonstrates, early- 
type systems typically have one bright PN (i.e., within 
0.5 mag of M*) for every ~ 5 x 10* of bolometric 
light, with redder galaxies trending towards lower val- 
ues of Q, and bluer systems approaching the theoretical 
limit mentioned above. But remarkably, recent mea- 
surements of ao.5 in the disks of late- type spirals have 
produced values of a in this same range, bet ween 1 and 



3 X 10" 



PNe L-^ 



Herrmann et all 120081 ). At first 



glance, the agreement between these numbers seems 
extraordinary. 

There are several factors which must be considered 
before interpreting this result. The first, of course, is 
that much of a spiral galaxy's luminosity comes from 
an extremely young stellar population that does not 
produce planetary nebulae. Specifically, if a galaxy has 
been producing stars at constant ra te over a Hubble 
time, then, according to the models o f iBruzual fc CharlotI 
(|2003l ). roughly half of its bolometric light should come 
from populations formed less than ~ 10* years ago. 
The luminosity from these young stars should dilute 
the galaxy's PN-to-light ratio, and cause our measure- 
ment of a to be lower than what would otherwise be 
predicted. 



rived from the exponential law. Observations of the 
brightest planetaries in late-type systems will there- 
fore overestimate the total number of PNe present in 
the galaxy. Fortuitously, this error is of the same order 
of magnitude as the luminosity dilution and works in 
the opposite direction. Thus, the observe values of a 
should probably not deviate very strongly from their 
nominal value. 

Interpreting the values of a seen in elliptical galax- 
ies is harder. PNe at the bright-end cutoff of the 
[O III] PNLF emit ~ 6OOL0 monochromatically at 
5007 A, and according to both models and observa- 
tions, no more than ~ 10% of the central star's flux can 
come out in this line (e.g.. lJacobvlll989l ; iMarigo et al.l 
|2004 | : ISchonberner et al.ll2007l ). Thus, the central stars 
bright PNe must be more luminous than 
in fact, in M31's bulge, 3 out of the 12 



of [O III 

~ 6OOOL0^ 

PNe analyzed bv I Jacobv fc CiarduUd l| 19991 ) have cen 
tral star luminosities in excess of 14, OOOLq . This is a 
problem, since in order to generate ~ 6OOOL0, a cen- 
tral star has to have a core mass of at least O.6M0 
(jVassiliadis & Wood 1994). Cores as m assive as this 
requi re high-mass (young) progenitors (jKalirai et al.l 
l2008i ). and such objects do not exist in elliptical and 
lenticular systems. 

There are two ways around this problem. The first 
is to hypothesize that objects at the bright-end of the 
[O III] PNLF are not PNe at all , but rather some 
sort of symbiotic star l|Soker|[2006D . Such a solution 
solves the luminosity problem at the expense of con- 
tinuity: why should such objects have the same max- 
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imum luminosity and same luminosity function nor- 
malization as normal planetary nebulae? The alter- 
native is to invoke a new channel of stellar evolution 
involving binary stars. Binary interactions on the red 
giant or asymptotic giant branches are unlikely to be 
the answer, since these encounters tend to truncate 
the stellar ev olution process befo re the creation of a 
massive core (|Iben fc Livid I1993D . However, encoun- 
ters on the main sequence, either through conserva- 
tive mass-transfer between primordial close binaries 
((McCrca 1964) or through the dynamical evolution 
of triple systems tPerets fc Fabrvck v '2009*) can cre- 
ate fie ld blue straggler stars, and, as CiarduUo et al.l 
(j2005l l have pointed out, blue stragglers have the proper 
number, core properties, and evolutionary timescales 
to evolve into [O III] -bright p lanet ar i es. In fact, if this 
scenario is correct, and the iKozail (|l962l ) cycle-tidal 
friction process is responsible for creating blue strag- 
gler stars, then many of the brightes t PNe we observe 
should have wide-binary companions (jPerets fc Fabrvckvl 
120091 '). 

How important is the contribution of blue strag- 
gler descendents to the PNLFs of elliptical galaxies? 
This is difficult to calculate, due to uncertainties in 
the initial mass-final mass relation at the low-mass 
end. However, if we assume that old stellar popula- 
tions are capable of producing ~ 0.55Mq cores, and 
that these cores can, at best, convert ionizing radiation 
into [O III] A5007 flux with 10% efficiency, then PNe 
from single stars can, theoretically, reach M5007 ~ —3. 
If the [O III] PNLF extends ~ 8 mag, then this implies 
that the brightest ~ 5% of planetary nebulae must 
come from an alternative channel of evolution. 



6 Prospects and Difficulties 

Over the years, attempts to model the PNLF have 
grown more and more sophisticated. Simple static 
models based on arbit rary distributions of central star 
masses (I Jacobvl 119891 have been replaced by simula- 
tions whMi^jasereaJistic galactic star-formation histo- 
ries l|Marigo et al.l200^ and hvd rodvnamic star-nebular 
interactions (jMendez et al.ll2008l ). These analyses have 
had some success. By including t h e effe cts of metallic- 
ity on the nebula, iDopita et al.l l| 19921 ) were able to 
reproduce the fadi ng of the PN LF cutoff in metal- 
poor populations l|Ciardullo et al. 200^ ) . Similarly , 
the hydrodynamical models of Mendez et al.l (|2008l ) 
have been able to accurately mimic the shape of the 
PNLF in young stellar systems, as well as the distri- 
bution of PNe in [O III] A5007-Ha-|-[N II] emission line 
space. 

The trouble comes when tryi ng to model t he be - 
havi or of the PNLF wit h age. iDopita et al] (|1992| ') 
and iMendez et al.l (j2008l ) both limited their analyses 
to the actively star-forming populations of the Mag- 
ellanic Clouds, with good reason. When considering 
single-star evolution, no reasonable initial mass-final 
mass relation can create the high mass cores needed 
to power M* plane t aries. This is illustr ated vividly by 
the simulations of iMarigo et al.l (j2004l '): ff the PNLF 
were populated solely by the remnants of single-stars. 



M* would fade by ~ 3 mag in the first ~ 6 Gyr af- 
ter the cessation of star formation. Clearly in order to 
reproduce the observed PNLF of old stellar systems, 
another channel is needed to create PNe with relatively 
high mass (Af > O.GMq) cores. 

Another problem with current PNLF modeling ef- 
forts is that none include the effects of circumstellar 
extinction. Most PNe are surrounded by a circum- 
stellar envelope of neutral gas and dust, and, as Fig- 
ure [6] illustrates, the extinction associated with this 
dust is an important contributor to the PNLF's cut- 
off. In star-forming systems, a typical [O III] bright 
planetary will be e xtinguished by A5007 ~ 0.7 mag 
([Herrmann fc Ciard uUo 200i|), and reddening values 
two to three time s this num ber are not uncommon 
(e.g., Mcathcring ham fc Dopit a 1 991ai ib|). A similar 

range of extinctions h ave been derived for the PNe of 

M31's bulge and M32 (jRicher et al.ll999l : lJacobv fc Ciardullol 
119991 ') ■ thus demonstrating that even in old systems, 
the effect of dust on the PNLF cutoff cannot be ne- 
glected. Of course, because extinction estimates for 
statistically complete sets of extragalactic PNe are dif- 
ficult to obtain, fully de-reddened PNLFs are not yet 
available to confront the models. Fortunately, this will 
soon change with the publicati on of the results from 
the MASH survey of the LMC (|Reid fc Parkei]l2009t ). 

Finally, the PNLF models to date have all ignored 
the information contained in the PNLF's normaliza- 
tion. A full description of the PNLF must not only 
explain its shape, but also the lower values of a ob- 
served in large, metal-rich elliptical galaxies, and the 
similarity between the a values of spheroid and disk 
populations. This may turn out to be the most diffi- 
cult challenge of all, if, as currently suggested, there 
is more than one scenario for planetary nebula forma- 
tion. 



7 Conclusion 

Most of the data on extragalactic PNe comes from ob- 
servations in the [O III] emission line at 5007 A. But 
the [O III] PNLF is just a single one-dimensional pro- 
jection of a more complicated function that distributes 
objects throughout a multi-dimensional phase-space 
that includes, for instance. Ha, H/3, continuum emis- 
sion at 4.5, 5.8, and 8.0 /jm, and the emission-line 
strengths of nitrogen, sulfur, neon, argon, and helium. 
Figure |4] displays a different one-dimension projection 
onto an axis which plots Ha, while Figure [S] shows a 
pair of two-dimensional projections onto planes which 
include the ratio of [O III] A5007 to Ha. A fully suc- 
cessful model of the PNLF needs to not only reproduce 
these projections, but also make predictions about the 
distribution of PNe in other phase-space dimensions. 
Similarly, the large homogeneous surveys now being 
performed in the Galaxy and the Magellanic Clouds 
have the potential to populate this PN phase space 
diagram as never before. The resulting confrontation 
between theory and observations is sure to generate 
new insights into the planetary nebula phenomenon. 
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